Stationary transport in mesoscopic hybrid structures with contacts to 
superconducting and normal wires. A Green's function approach for multiterminal 

setups. 
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We generalize the representation of the real time Green's functions introduced by Langreth and 
Nordlander [Phys. Rev. B 43 2541 (1991)] and Meir and Wingreen [Phys. Rev. Lett. 68 2512 
(1992)] in stationary quantum transport in order to study problems with hybrid structures containing 
normal (N) and superconducting (S) pieces. We illustrate the treatment in a S-N junction under 
a stationary bias and investigate in detail the behavior of the equilibrium currents in a normal 
ring threaded by a magnetic flux with attached superconducting wires at equilibrium. We analyze 
the flux sensitivity of the Andreev states and we show that their response is equivalent to the one 
corresponding to the Cooper pairs with momentum g = in an isolated superconducting ring. 
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I. INTRODUCTION. 

The superconductivity and its implications is among 
the most interesting phenomena in the realm of con- 
densed matter physics. While the microscopic mech- 
anism leading to the pairing instability in the high-Tc 
materials remains not yet fully understood, the general 
framework provided by the BCS theory^ consistently ac- 
counts for superconductivity in normal metals. Remark- 
ably, this seems to be even true in the context of low 
dimensional systems of mesoscopic scale. ^'^ 

BCS theory provided the basis of the seminal paper by 
Blonder, Tinkham and Klapwijk (BTK)^. In that work, 
the stationary transport properties of a superconductor- 
normal metal (S-N) junction and the subtle mechanism 
of the Andreev reflection leading to the effective Cooper 
pair tunneling through the junction was first analyzed. 
A similar description was followed in the study of S-N-S 
structures^'i^iii^ii^ and later formulated in terms of multi- 
channel scattering matrix theory in Ref [lo|. BCS theory 
has been also the basis for the study of stationary trans- 
port in unbiased S-N-S junctions due to the Josephson 
effecljiiiiii^ii^iii^ as well as the AC Josephson effect un- 
der bia9i>Ha5i6a7_ 

The non-equilibrium Green's function formalism^- is 
a powerful technique to study quantum transport in co- 
herent regimes. In the context of microscopic models for 
mesoscopic structures it was first introduced by Caroli et 
al^ and later elaborated by other authors . ^^i^^'^^i^^'^*^!^^ 
That approach was also represented in the Nambu 
formalism to treat S-N and S-N-S junctions.— i^^i^i^ 
The formal equivalence between non-equilibrium Green's 
function and the scattering matrix formalism to the 
quantum transport has been analyzed for the case of 
normal systems without many-body interactions under 
stationar}*22, and time-periodic drivingi^^ 

The representation of the non-equilibrium Green's 



functions introduced by Langreth and Nordlander— 
is particularly useful to derive compact equations for 
the currents along the different pieces of a mesoscopic 
structure.—"^ In the present work, we employ that rep- 
resentation in the case of hybrid multiterminal structures 
containing superconducting elements that are modeled 
by BCS Hamiltonians. 

Instead of working in Nambu's space, we derive a cou- 
pled set of Dyson's equations for the normal G^'^ (w) and 
Gorkov's F^''^{u!) retarded (i?) and lesser (<) Green's 
functions. As in Refs .^^'^'^i^^ , we "integrate- out" the de- 
grees of freedom of the external wires (reservoirs) and, 

by introducing auxiliary hole propagators g ^(i-ij), we 
reduce the problem to solving the Dyson's equation for 
the usual normal Green's function with an effective self- 
energy. As in Refsi^i^^i^i, the latter describes the scat- 
tering events due to the escape to the leads, but in the 
present case, it contains a component related to the mul- 
tiscattering processes involved in the Andreev reflection. 
The final expressions for the currents have a c omp act 
structure that formally resemble those of Ref. l23| for 
normal systems. 

Sections II and III are devoted to explain the theoret- 
ical treatment. We derive expressions for the currents 
and we show that the transmission function of a biased 
system contains a normal plus an Andreev contribution. 
In Section IV we illustrate the approach in the simple 
well known case of a two terminal setup with a linear 
system in contact to one normal and one superconduct- 
ing wires under bias and we show its equivalence with 
BTK description. In Section V we employ the formalism 
to the study of the behavior of the equilibrium currents 
of a normal metallic ring threaded by a static magnetic 
field with several attached normal and/or superconduct- 
ing wires. We address several interesting physical ques- 
tions like the minimal conditions for the development of 
Andreev states within the superconducting gap, the flux 
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FIG. 1; (Color online) Sketch of the setup. The central grid 
represents the central finite system. The area enclosed by 
this system is threaded by a static magnetic flux $. The 
N- and S wires are, respectively, indicated with open and 
filled lines. The arrows represent the contacts between the 
different systems. In each case, the parameters of the ensuing 
Hamiltonians are indicated. 



sensitivity of these states and the possibility of anoma- 
lous flux quantization induced as a consequence of the 
proximity effect. Section VI is devoted to summary and 
discussion. Some technical details are presented in the 
appendices. 



with a =1, i, and being = for the N-wires. The size 
of the wires approaches the thermodynamic limit {La — > 
oo), i.e., the wires act as macroscopic reservoirs, with well 
defined chemical potential and temperature. We model 
the central system by a tight-binding Hamiltonian in a 
finite lattice of L sites with nearest-neighbor hopping. 
We consider the possibility of a static magnetic fiux <i> 
threading this system, which introduces a dependence 
on $ in the hopping matrix elements: 

L 

+ E s^c^Q,., (3) 

l=l,(T 

where (W) denotes nearest-neighbor sites. The Hamilto- 
nians for the contacts read: 

Hca = ~^caE(4<:„,o-Cicc,o- + H.C.), (4) 

CT 

which describe hopping processes between the sites jca of 
the wires and the sites lea of the central system at which 
the wires are attached. 



II. THEORETICAL TREATMENT. 
A. Model 

We introduce microscopic models for the different 
pieces of the setup, which consists in a finite normal sys- 
tem of non-interacting electrons in contact to M infinite 
superconducting (S) or normal (N) metallic wires (see 
Fig. 1). The full system is described by the following 
Hamiltonian: 

M 

H = Hcen + E^^" ^ Hea)^ (1) 
a=l 

where Ha denote the Hamiltonians of the wires, while 
Hea the corresponding contacts establishing the connec- 
tions between these systems and the central one. Al- 
though long-range superconducting order does not take 
place in strictly onc-dimcnsion (ID), for simplicity, we 
consider ID tight-binding BCS Hamiltonians with local 
s-wave pairing for the wires. This is a rather standard 
assumption (see Refs. iiiHilUllliillll) and the 
general treatment can be easily extended to multichannel 
wires and more general symmetries of the superconduct- 
ing gap. Concretely: 

Lc 

Ha = ~Wa 2 (cI,'tCj= + 1,<t +ii'.C.) - (2) 

ja=i,<y 



B. Currents. 

The electronic current, in units of e/h, flowing through 
a given bond {I, I') of the central system is: 

^^Re[w,,i{<P)G<,^A^)], (5) 
while the current flowing through a given contact is 

/OC J 
^^Re[WeaG<_^_A^)], (6) 

being 

Gfj,Jt,t') = ^{cl{t)c,,it')}, (7) 

and Gfi, ^(10) the corresponding Fourier transform in t — 
t'. 



C. Evaluation of the Green's functions. 

In previous literature, the evaluation of the Green's 
functions for hybrid structures described in terms of 
tight-binding and BCS Hamiltonians has been carried 
out in the framework of the Nambu formalis m^^'^^1^^1^^ . 
We briefly present bellow an alternative and equivalent 
representation, which will allow us to analyze from a dif- 
ferent perspective the physical processes involved in the 
phenomena of Andreev reflection and the development of 
Andreev states within the superconducting gap. 
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We define retarded normal and Gor'kov Green's func- 
tions: 

G«,,,(t,t') = -^e{t-t'){{c,At),cl,Jt')}), 

where {., .} denotes the anticommutator of the corre- 
sponding operators and t =i) i =T- 

It can be verified that the equations of motion for these 
functions are coupled and read: 

j" 
j" 

The spacial indexes extend over the coordinates of 
the whole system. For coordinates on the wires 

For coordinates on the central system: Sjj' = 

for {1,1'), being nearest neigh- 
bors within the L-site lattice, ejj' = ^'^'^'■j'^ij' 
and Aj = 0. For coordinates on the contacts: ej,j> = 
-Wca{5j.i^^5j' ^-j^^ + Sj.j^^Sf ^i^^) and A; = 0. 

As usual, it is convenient to eliminate the degrees 
of freedom of the wires. Such a procedure defines 
self-energies for the Green's functions with coordinates 
belonging to what we have defined as the central 
system, We summarize it in Appendix |A] for the 
present problem. The result is that the retarded Green's 
functions with coordinates on the central system can be 
expressed as elements oi L x L matrices and the ensuing 
Dyson's equations read: 

[f{u:)r^F^{u:) + E/«^«(c.)G«(c.) = 6, (8) 

where T.'^'^i, {uj) = 5i,vY.Ju..^7''"i'-^). with = 
g, f. The explicit evaluation of these functions is summa- 
rized in Appendix B. The have introduced the retarded 
Green's functions g^{Lo) and g [lo), whose corresponding 
inverses are: 

[f{uo)]-' =Lui + £(-$) - tff^^'iLu), (9) 

where £($) contains the matrix elements of the Hamil- 
tonian Hcen- In the case that all the wires are nor- 
mal (Aq = 0,Va), the function g^{Lo) is the exact re- 
tarded normal Green's function of the coupled central 
system, while Y.ff^^{uo) = thus = 

[g^{—uj)]*, which indicates that 'g (cu) is a propagator 
related to the dynamics of the holes. 
The second equation ([5]) can be casted: 

F^icj) = -fiu;)±f^-"ic.)G^{c.). (10) 



Substituting ((TO)) in the first equation ([8]) the formal so- 
lution for the normal Green's is obtained: 

[G^H]-i^..i-e-Sf^/H, (11) 

where we have defined an effective normal self-energy: 

tfffiuj) = SS»^^(t^) + E9-^'^(c^)t'^(w)E^s^^(t^). (12) 
The lesser counterpart of pT|) is, thus, written as: 

G<(u;) = G^(c.)E</.;)G^(^), (13) 

being the advanced Green's function G^{uj) = [G^{uj)]^ 
Using Langreth rules^: {BC)< = B^G<+B<G^ in the 
definition of (|12p . it can be shown that 

x[f{uj)tfs^^{u:)^f \uj)tfs^<{u:)]. (14) 
Using the lesser counterpart of ([9]): 

t[u^)^f{u)tff'<{u)f{u), (15) 

the lesser effective self-energy S^yy(tj) can be fully 
expressed in terms of the bare ones, S^"^ '^(^) = 
ifa{Lo)T^'^ (w), with v,v' = g,f, which depend on the 
temperature Tq, of the reservoirs through the Fermi func- 
tion /a(w): 

+ Y.klA^)Y.fJr<{u^)Kt,^p{u;), (16) 

ol' 

with K^pioj) = Ss/^«(a;)5«^^;^^(c.) and A^ Jc.) = 

[K^ piuj)]* . Alternatively, the above expressions can be 
also directly obtained after some algebra from the lesser 
counterpart of ([S]), as indicated in Appendix C. 
At equilibrium, it is satisfied: 

^^//,a,/3 = if{^)^eff,a,fi{uj), (17) 

being 

M 
a' = l 

which implies: 

G<,Ju;) = fiu;)[Gt,Ju;) - Gj',,Ju;)]. (19) 
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Before closing this section, let us emphasize the formal 
equivalence between Eqs. pT|) and and the represen- 
tation of Ref. [2ll[23H23 . In the present case, the effective 
self-energies (|12p and p^ . however, have a more com- 
plicated structure when the leads are superconducting. 
In particular, they contain the normal terms Y,^^'^'"^ (lu) 
that represent the normal "escape to the leads" of sin- 
gle electrons, as well as terms involving multiple scatter- 

ing processes, mediated by the hole propagators g ' (uj). 
The latter act not only locally, but also extend along the 
different positions of the sample that are in contact to 
superconducting wires. 



III. STATIONARY CURRENTS AND 
TRANSMISSION FUNCTIONS. 

Being able to evaluate the lesser Green's functions, we 
are now in the position to evaluate the currents (O and 
We recall that a biased setup with several super- 
conducting wires defines, in general, a time-dependent 
problem.— In this work we are interested in the sta- 
tionary transport. Thus, in what follows we shall derive 
expressions for the currents in two situations: (i) A bi- 
ased setup with a voltage difference between the S- and 
the N- wires, being all the S wires at the same chemical 
potential. In this case, currents flow through the con- 
tacts as well as along the central system, (ii) The second 
situation corresponds to all the wires at the same chemi- 
cal potential, in which case, there are no currents flowing 
through the contacts and there exists only the possibility 
of equilibrium currents along the central structure when 
it is threaded by a finite magnetic flux. We present below 
general exact expressions for the currents and we shall 
address separately the two different cases in the next two 
sections. 

Using Dyson's equation for the lesser Green's function, 
the expressions ^ and (O cast: 

Ji^i, = -2 V / -^Re[u;i,,i($) x 

G« (^)S<//,„^ {uj)Gl^ {uj)i (20) 
for the current along a given bond (Z, V) and 

M f + oo 1 

+ S«/,„^(C.)G<^_,^^^^(^)], (21) 

for the current along the contact to the wire a. Details 
for the derivation of the latter equation from ([6|) follow 
the same lines as in Refs. [23ll2i (see e.g. Eq. (5) of Ref. 
[23h . using the normal Green's functions ITTI) and ([13]). 



A. Equilibrium currents. 

When the central system is attached to wires at 
the same chemical potential fi, there is no charge flow 
through the contacts to the reservoirs. Nevertheless, if 
the central system is threaded by a finite magnetic flux, 
equilibrium currents can flow within this system. For a 
given bond {1,1'), the equilibrium current reads: 

TtJ,{u;) = -2Re{w,,imGti'A^)-Gf^,,Acj)]} 

M 

= 2Im[ ^ re//,„/3(tj)«;r,i(*) 

xGj^i^^Jij)Gt^^,Jij)], (22) 

where we have used the equilibrium identities (jl7[) and 
(HH), while re//,a/3(w) is defined in Eq. For $ = 0, 

the result T^l, {u!)\,s>=o = is obtained by noticing that 
the function within [. . .] of the above expression is just 
the real function — 2Im[w;/^;(0)Gf-j,(a;)|$=o]- 

B. Non-equilibrium currents. 

We consider Ms S-wires at /io, = /i and M_\f = M — Ms 
N-wires with a voltage difference V with respect to 
the superconducting ones. Following Ref. ^15" we take 
/ict = /X in the Hamiltonians Ha for the N-wires and en- 
close the bias V in the corresponding Fermi functions. 
We also consider that all the wires are at the same 
temperature. Therefore, for the N-wires: 'E,^^''^{uj) = 
if{Lj-V)ra{uj) and T:fJ'<{uj) = if{uj + V)Ta{uj), where 
Taioj) = rs3(w)|A„=o, while for the superconducting 
ones: T.''/'<{lu) = i/(tj)r^'''(w), with = g,f. In 
order to derive the expressions for the currents it is use- 
ful to express the effective lesser self-energy as follows: 

^fff.aA^) ^if{^^)^eff,a,pH+ (23) 

t[fiu;-V)-f{iu)]Sa,f^ J2 Sa,a'rmu;) + 

where Tef f^cpi^jj) has been defined in Eq. pS]) . 

The final expression for the non-equilibrium current 
along a given bond of nearest neighbors {I, I') is: 

^[f{u;-V)-fiu;mcj). (24) 
-oo 27r 

In the case that, in addition to the bias V, the central 
system is threaded by a magnetic flux, we should add 
to the previous expression the equilibrium contribution 
Ji^i defined in the previous subsection. J^l, is due to the 
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internal currents of the single-electron orbits of the finite 
system that are twisted by the static flux. Instead, the 
origin of the non-equilibrium contribution is a net parti- 
cle flow between reservoirs through the central structure. 
For this reason, the non-equilibrium component depends 
only on the spectral properties within the energy window 
[/i,/i-f y], while the equilibrium one formally depends on 
the spectral weight of all the quantum states bellow fi. 
The transmission function contains two contributions: 

Tu,{u)=TI:A^)-T[:A~lo). (25) 

The first one is the normal transmission function: 

Mm 

cr,aeAf=l 

ImK,,('i>)Gfi,„,.(^)G^^„,iv(a;)], (26) 
and the second one is the Andreev transmission function. 

Mat 

Tr,,iio) =-2 J2 r^^'Hx 

cr,QGAA=l 

ImK,K*)AL,.('^)Ai/',.('^)], (27) 

where the a ^ J\f denotes summation over the nor- 

mal wires, while A, „^^(tj) = E/3 G^^^ ,^A|^„(cj) and 
fi 

A-Q,/'(w) = [A;,_^(w)]*. While the normal transmission 
function depends on the rate at which electrons can be 
emitted at the normal reservoirs r^^(a;), the Andreev 
transmission function depends on the rate of emission of 
holes (we recall that r£-''(cj) = rs9(-w)). The Andreev 
component depends on the multiple scattering propaga- 
tors A; ^{lo). Instead, the normal component depends 
on the usual ones G^i^^ (t(^)- For ^ vanishing supercon- 
ducting gap, Tj" , {ld) = 0, and only the normal compo- 
nent survives. 

Analogously, the currents through the contacts can be 
written as: 

/ + OC 7 
^[f[u:-V)- f{uo)]T^[i^), (28) 
-oo 2tT 

with the transmission function also containing two com- 
ponents: 

T^iu) = T2{lo)-T-{~u). (29) 
The normal transmission function reads; 

M 

a, (3=1 
Mm 

Y C(^)Im[Sf^(a;) 

a'eAf=l 

xGf^^i^^,^A^)Gf^^,i^^jLo)]}, (30) 



while the Andreev transmission function is: 

M Mm 

xIm[Al„,Ju;)A^,^^Ju;)Gt^,^^Ju;) + 
Sf//.a0(^)A^„,.,,(^)Af,,^^,,(c.)]. (31) 

IV. A LINEAR BIASED SETUP WITH A 
SINGLE SUPERCONDUCTING WIRE AND A 
SINGLE NORMAL WIRE. 

In this section, we shall explicitly write down the pre- 
vious expressions for the case of a setup with two wires: 
one superconducting and the other one normal, which we 
denote, respectively, a — N and a = S. This will allow 
us to show that we are able to recover BTK's description^ 
for the transmission functions of a simple tunneling junc- 
tion. 

In this case: i;f^y^^^(a;) = SaA^a.N^^^'^i'^) + 

Sa,s^s^''^'{^)gil isi^)'^s^'''^{^)- The total transmission 
function evaluated at the contact with the A^-wire is 
Tis!{Ld) = T'^{lo) — T^{~llj). The normal component is 
given by Eq. (|30p . which in this simple case reduces to: 

TU^)=J2^%^(.^)\Gt,isA^)\''^'eff.si''y (32) 

Notice that we recover the well known structure for the 
normal transmission function in terms of Green's func- 
tions originally pointed out by Fisher and Lec^". In 
the present case, the function ^Ijf gi^) = ^s^i'^) ^ 
2Im[E|/'^([j)^;^ i^{Lu)I]g^'''^'{uj)] contains the usual term 
r|r'(cj), which depends on the normal density of states of 
the superconducting lead, as well as a multiple-scattering 
term that depends on the hole propagator 'gf^ (uj) and 

the anomalous self-energy of the wire Ylg^'^{u!). The An- 
dreev transmission function reads: 

T^iu;) = -J2^fJiu;)\A'^^^Ju;)fT^j^iu;), (33) 

which actually has the formal structure of a reflection 
process represented in terms of Green's functions. Fur- 
thermore, it depends on the emission rate for holes in the 
normal wire rjj{uj) = F^ {—^) and it contains a multiple 
scattering kernel: 

A^,^,,(a;) = Gf„,,,^,(.;)A«^(a;), 
A§^^iu;) = ^■(c.)g«,„(.;). (34) 

After some algebra, it can be verified that: T^{lu) = 
-T^(cj) = T^iuj) and TJ^^iu;) = -Tg{Lj) = T"(tj), in 
consistency with the continuity of the current. 

In order to benchmark the above representation, we 
present results for the central system being a linear one- 
dimensional junction with a barrier of height Eq as in 
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FIG. 2; (Color online) Benchmark against BTK theory. 
Transmission functions T"{uj) (dashed black lines) and T"(lj) 
(red solid lines) in the lower panels and the total transmis- 
sion T{uj) — T"{uj) — r"(tj) in the upper panels for a junction 
described by the Hamiltonian (|35|l . Left and right panels 
correspond to Eo = 0, 1, respectively. Other parameters are 
wn ~ ws — w = 1, jj, = and As — 0.2. 

BTK's paperi. (see also Ref. i): 



Hcen = y^[c[^Q + i,g + H.c] 

1 

+ E (35) 
i=-i 

with III — "^rj cj ^ci,cr aiid = — /i + Si^qEq- For such 
a system, it is easy to verify that the expressions for 
the transmission functions corresponding to a given bond 
{I, I + 1) are: 

(T 

a 

x\Ks,nA^)?'^n{^)- (36) 

It can be proved that this functions also satisfy 
r/;;^^(tj) = T"(a;) and r°j+i(a;) = T"(u;), in agreement 
with the conservation of the current. 

Numerical results for the functions T"(a;) and T°'{lu) 
are shown in the lower panels of Fig. [5] The correspond- 
ing total transmission T{uj) is also shown in the upper 
panels for Eq = and Eq = 1. The picture presented in 
BTK's paperi is identified through T"(w) ^ 1 - B{E) 
and T^ito) -A{E), with A{E),B{E) defined in Ref. 
|3. The lower panels of Fig. 2 should be compared with 
Fig. 5 of Ref. |j. It is worth noticing, in particular, the 
fact that T°'{uj) is sizable within the gap, while in the ab- 
sence of a barrier {Eq = 0), T'^(w) -> -1. Thus T(w) - 2 
for < a, (see upper panels of Fig. 2 and compare with 
Fig. 7 of Refi). 
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FIG. 3: (Color online) Sketch of the setup. The central system 
is a ring threaded by a magnetic flux in contact to supercon- 
ducting and normal reservoirs at the same chemical potential 
/i. The only non-vanishing current is the equilibrium current 
along the circumference of the ring. 



V. FLUX SENSITIVITY OF THE 
EQUILIBRIUM CURRENTS IN A RING. 

We now turn to the setup without bias voltage {V 
We consider the simple case sketched in Fig. [31 where 
the central system corresponds to a one-dimensional ring 
threaded by a magnetic flux <I>, i.e. Hcen = Hrmg, being: 

L 

Hring = -W ^ (e"'*/^c|^^Q+i,^ + iJ.C.) 

1=1,(7 

L 

l=l,cr 

where $ is expressed in units of 27r$0j being $0 — e/h 
the elementary quantum. We take the lattice constant 
a = 1 and we impose the periodic boundary condition 
L-f 1 = 1. 

An isolated normal ring under a magnetic flux, sup- 
ports a persistent current with a periodicity equal to $0, 
as a consequence of the sensitivity of its energy levels 
with the threading flux. When normal metallic wires are 
attached to the ring, inelastic scattering effects are intro- 
duced which decrease the magnitude of this equilibrium 
current. However, its qualitative behavior, in particular, 
the periodicity with the flux is expected to be the same 
as in the case of the isolated ring, provided that the in- 
elastic scattering length introduced by the coupling to 
the external wires satisfies > La. For < La, this 
current is, instead, expected to vanish. This is because, 
for a short enough ring such that > La, the effect of 
the coupling to the wires is essentially the introduce tion 
of a finite lifetime in the energy levels, without affecting 
their fiux sensitivity. 
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In the case of an isolated superconducting ring with s- 
wave pairing, Byers and Yang2^ have shown that the pe- 
riodicity of the flux-induced persistent currents is $o/2- 
This is again a consequence of the sensitivity of the en- 
ergy levels, this time combined with the fact that the 
structure of the wave function corresponds to an ensem- 
ble of Cooper pairs, instead of one of single electrons. Hy- 
brid isolated S-N piecewise rings have been also studied 
and the conclusion is that the periodicity of the persis- 
tent currents experiences a crossover between $o/2 and 
$0) as the length of the superconducting piece becomes 
shorter than the superconducting coherence length ^t—i^. 

On the other hand, a conductor between two supercon- 
ductors forming a S-N-S structure is known to support 
Andreev states within the superconducting gap. In par- 
ticular, such states are expected to develop for a ring 
with attached superconducting wires and it is interesting 
to study the flux sensitivity of these states, which should 
define the behavior of the equilibrium currents. It is also 
interesting to investigate which is the minimum number 
of S-wires needed to develop Andreev states. Further- 
more, recent studies suggest that the vortex excitations 
of a superconducting state can exist within a normal con- 
ductor sandwiched between two superconductors^ due 
to the proximity effect. It is, therefore interesting to in- 
vestigate whether it is possible that proximity effect in- 
duces also a flux periodicity of $0/2 in a normal ring due 
to the attachment to S-wires. 

In order to address these issues we analyze the behavior 
of the function T'^'^iuj). Because of the continuity of the 
charge, this function is independent of the bond 1,1 + 1 
along the ring chosen for the evaluation of Eq. ([22|) . 
Thus, the latter expression can also be written as follows: 



2w 



-ii>/L 



1=1 (J,a,p 



x[Gf,+i.,H-[Gf+i,,,H]*]}. (38) 



In what follows, we analyze different configurations of 



A. Each site of the ring in contact with a wire. 

Let us first consider the simple case of a ring in con- 
tact to wires in a configuration that does not break the 
periodic translational invariance along the circumference 
of the ring. Such a configuration corresponds to L iden- 
tical wires (N or S), each one in contact to a single site 
of the ring. The retarded Green's function can be easily 
evaluated in this case. The result is: 



771 — 



(39) 



with km = — TT + 2m7r/L, m = 0, . . . , L — 1, and £m{^) = 
~2w cos{km + ^ I L) , where, for simplicity, we have taken 
fi = 0. The effective self-energy is: 



H = S]^^'^(c.) - j:^f'^iu;)g^{co)J:f^'^ico), (40) 



where the second term vanishes for N-wires. The hole 
propagator of this term is: 



+ £„(-$) -E//'«(C^)' 



(41) 



Transforming the right hand side of psp to the reciprocal 
space, it reduces to: 



T^'^iu;) = I Xi «™('i>){-2Im[G^,(^)]}, (42) 

m=0 

with Wm(<i>) — 2wsm{km — ^ / L) ~ dem{^)/dk,n being 
the velocity corresponding to the m-th energy level. 

In the limit where the coupling to the wires vanishes, 
the above expression reduces to the transmission function 
of an isolated ring: 



m— 



m)- (43) 



For N-wires or for S-wires and energies such that 1^1 > 
A, a similar expression is obtained: 



z;,„(<i>)Im[E^{/'«(c.)] 



L 



m=0 



(44) 

where the 0-function applies only for the case of a S- 
wire. The above expression corresponds to a sequence 
of Lorenzian functions centered at energies ~ £„(<&) + 
Re (£,„($))] with width ~ Im[EJ^//'«(e„($))]. 
The latter parameter defines the lifetime of the levels 
of the ring due to the coupling to the reservoirs. 

The periodicity of these currents as functions of the 
flux is $0, which corresponds to a shift $/L = 27r/L, that 
is equivalent to a relabeling of the reciprocal points km- 
For S-wires and |w| < A, the functions V^'" [u) = 0, thus 
lTsi\T,'^f'^{Lo)\ = 0, and the only spectral contribution to 
T'^'^{bj) is due to the eventual development of Andreev 
states. The energies of these states is determined from 
the poles of the function which implies finding 

the roots of the function: 

\{u;) = CO -£„,($) -Re[I]9S'«(a;)] 

-Re[Y.^f^"{oj)i:f^^\u:)]Re[g^{oj% (45) 



where 



5^(c.) = e(A-|c.|)- 



1 



UJ + £-„($) + W 
with £,„($) ~ Emm + RepfS (£„($))]. 



(46) 
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Approximating Re[E''''' (w)] ~ Rep^''' (£±m (*))], the 
solution casts the following roots: 

E i{<t>) ^ e-($) ± 

^[e+(<i>)]2+Re[E9/(£„($))S]/9(e,„($))], (47) 



with 



+ _ £r«W±£-m(^) 



(48) 



while the corresponding quasiparticle weights are: 

— TT — TT 



7± - 



Replacing in (|42|) . the final result for the transmission 
function within the superconducting gap is: 



27re(A- \lu\) 



L-l 



s— it ,m 



=0 



For \uj\ < A: 



Re[ESS(tj)] = W7(w), 



being 



7(w) 



2wl 



4wi 



A2 



(50) 



(51) 



(52) 



Therefore: 



- /?[£,„($) -£-,n($)]± 
V/32[e™($)+e-™($)]2+72A2, (53) 

being /3 = (1 + 7(e™($)))/2, and 7 - 7(em($))- 

Remarkably, the expression (j50p with the energy given 
by l|53p coincides with the expression for the persistent 
currents of an isolated ID BCS tight-binding ring with 
hopping 2/3u;, gap 27 A and pairs with total momentum 
q = (see Ref. ISOl ). In other words, the flux sensitiv- 
ity of the Andreev states in our problem is exactly the 
same as that observed in an isolated BCS ID ring with 
pairs of momentum q = 0. The fact that only pairs with 
momentum 5 = contribute implies that the periodic- 
ity of these currents is just the normal periodicity of a 
flux quantum $o- These currents do not show the $o/2 
periodicity, typical of a true superconducting ring, since 
the origin of that behavior is a change in 27r/L of the 
total momentum q of the Cooper pairs. The renormal- 
ization factor (3 for the hopping parameter within the 
ring, which determines the velocity and, thus, the 
amplitude of the currents, depends on the superconduct- 
ing coherence length of the wires, ~ A/2w, as well as 
on the tunneling ratio through the contacts, controlled 
by the parameter Wc- Its magnitude is large for energies 
close to the edge of the gap |em((/))| ^ A. 



B. A single S-wire attached to the ring. 

Let us now consider a single superconducting wire at- 
tached to the ring. 

As before, we must consider separately the contribu- 
tion from states with energies within and away from the 
superconducting gap. To analyze the spectrum for en- 
ergies \u}\ > A, it is convenient to write the retarded 
Green's function as follows: 



(54) 



(49) being 



m=0 



9lM = 



1 



(55) 



andEf^.^^Ju;) = J:^^{u;)+^^J {u;)gl^^i^Ju;)j:fj{cj), with: 



9^,1' A) 



1 

m=0 



-ifc„(i-i')-0 



9UA, 



(56) 



Substituting in (j38p . the transmission function reads: 
T^.(^) ^ ?^M^^\„($)^„(^), (57) 

m=0 



being 



^m(w) = 



re//,a,a(t^)|gg^ (t^)l^ 
\^~^fffJA9LlJA\ 



2 ' 



(58) 



which results in a Lorcntzian-type proflle as in the case 
ofEq. (gl. 

As in the case considered in the previous subsection, 
for < A, lm[Ts'^''''^{uj)] = 0, and Andreev states can 
develop within the gap. In order to determine the en- 
ergies of these levels, it is convenient to consider the re- 
tarded Green's functions gf^i,{uj) and gj^i,{uj), deflned in 
Eqs. ([9]), which in the present case are the solutions of 
the following Dyson's equations: 

g« (c) = glAA+9tM^T''A)9LMiA, 

gf,,{A - fiAA+9LMK^'''iA9l.AA-m 

Within the gap, these functions have, respectively, quasi- 
particle and quasihole states, behaving as follows: 



9f:vA) 



0(k|-A) 



X - e 



LO- Em ($) + iV 



L-l 



L 



qU}+ e_.m ($) -t-i?7 
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being £„, ($) - e™($)+CRe[Sfl;i?(£„($))]/L, where C = 
2 for $ = i^TT with K integer while C — 1 otherwise, and 
Z,„ = -7r{|l-C9Re[S99,fl(^)]/a,^|^^^^^/L}-i. in what 

foUows, we shall approximate ^ — tt, which becomes 
exact in the limit L — > oo. 

The full retarded Green's function is, in turn, deter- 
mined from: 

X5L,.J^)S^^'^(^).9/l,r(^)- (61) 

As in the previous section, the ensuing solution has a 
quasiparticle BCS-like structure: 



L 



u - E?„ ($) + iri 



b(62) 



with E^{^) given in with 7 cx 1/L and given 

in Eq. (01]). 

Therefore, for a single superconducting wire connected 
to a large enough ring, Andreev levels tend to coin- 
cide with free particle and hole energies: and 
— e_m(<i>), respectively, provided that |em(<i>)| < A, 
k-m('i')| < A. The corresponding transmission function 
is formally given by Eq. ([50|) . 

In conclusion, a single superconducting wire attached 
to the ring generates the same qualitative behavior as L 
superconducting wires attached in a translational sym- 
metrical way, but the effect is 0{1/L) and tends to be 
negligible as i — > 00. 



VI. SUMMARY AND CONCLUSIONS. 

We have presented a representation of Keldysh Green's 
functions for stationary transport problems in systems 
with superconducting and normal components. As most 
of the relevant observables, like the currents, depend on 
normal propagators, we have worked with Dyson's equa- 
tions in order to eliminate the anomalous ones. This 
procedure has been carried out by defining auxiliary hole 
propagators and effective self-energies that contain mul- 
tiscattering terms. In the resulting representation, the 
Green's functions exhibit the same structure as in nor- 
mal systems. This allows for the derivation of simple and 
compact expressions for the currents and the transmis- 
sion functions, that are similar to the ones presented in 
Refs. [23112^ for normal systems. 

We have presented general expressions for the currents 
in stationary conditions, distinguishing two situations: 
biased systems where transport is induced by a voltage 
difference and equilibrium currents induced by a static 
magnetic flux. In the case of biased systems, we have de- 
fined normal and Andreev transmission functions and we 
have compared them with results obtained in the frame- 
work of previous formalisms, in particular, the one pre- 
sented by Blonder, Tinkham and Klapwijk. 



We have, finally focused in the study of the behav- 
ior of the equilibrium currents in a tight-binding normal 
ring with attached superconducting wires. These cur- 
rents result as superpositions of the currents of all the 
states of the ring with energies e,n(<I') bellow the chem- 
ical potential of the wires, in which electrons circulate 
with velocities Vm = d£m{^) / dkm- 

Our main conclusions on the qualitative behavior of 
these currents are the following: (i) The states with en- 
ergies lying away from the energy window defined by the 
superconducting gap present an identical qualitative be- 
havior as those of rings attached to N wires. In partic- 
ular, they have a periodicity of $0 functions of the 
external flux. The spectral profile related to these cur- 
rents is a collection of Lorentzian functions which implies 
a decrease in the amplitude of the current due to inelastic 
scattering effects via the escape to the leads. 

(ii) The states with energies within the superconduct- 
ing gap of the wires, behave as isolated in the sense that 
the spectral weight related to them consists in a collec- 
tion of delta functions, indicating the lack of inelastic 
scattering effects. The positions of the energy levels is, 
however, affected by the proximity effect and they are 
organized in a structure that replicates the quasiparticle 
spectrum of a BCS tight-binding superconducting ring 
with Cooper pairs of momentum 5 = 0. The effective 
BCS tight-binding parameters are the hopping, which is 
the bare hopping of the ring renormalized by a factor (3 
and a gap, which is the gap of the superconducting wires 
renormalized by a factor 7. The renormalizing factors 
depend on the superconducting coherence length of the 
wires and the degree of coupling between the wires and 
the ring. The latter effect is controlled by the strength 
of the coupling between these systems as well as on the 
number of attached wires. For a single attached wire, 
it is 0{1/L) and, thus, not significant for large enough 
rings. 

(iii) Although the proximity effect induces Andreev 
levels that replicate the structure of quasiparticle states 
of a superconducting ring within the energy window de- 
fined by the superconducting gap of the wires, these 
states correspond only to the subspace with winding 
number q = Q. Since the periodicity in ^^0/2 of the per- 
sistent currents in superconducting rings is explained by 
a shift in the winding number q commensurate with the 
reciprocal lattice of the i-ing,^*'^° the restriction of the 
subspace with q — Q does not allow for such a mecha- 
nism. The consequence of this rigidity is that Andreev 
states have the same periodicity $0 as the states of the 
normal ring. Let us, however, mention that the rigidity 
of the winding number could be due to the rigid BCS 
mean field approximation considered to model the exter- 
nal wires. There exists the possibility that a more flexible 
model allowing for spacial fluctuations of the parameter 
A within a region of the external wires that is close to 
the contacts could also permit fluctuations in the winding 
number q of the induced Andreev sates within the ring. 
A possibility to explore this mechanism is by recourse to 
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a self-consistent approximation similar to that of Refs. 
[Til and [13. 
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APPENDIX A: ELIMINATING THE DEGREES 
OF FREEDOM OF THE RESERVOIRS. 

We summarize the procedure introduced in Ref. [2^ . 
[23l[23 to eliminate the degrees of freedom of the external 
wires in the Dyson's equation for the central system. 

It is convenient to change the basis in Ha as follows: 



N. 



a + . n 



(Al) 



ri=0 



with kn^a — n-K I [Na + 1), n = 0, . 



. , Na, which leads to: 



Hn 



n— cr 

No. 



(A2) 



being efe„^ = 



n=0 

-2WaC0Skn^a — Mai ^^'^ 



No, 



(A3) 



n=0 cr 



being Wa,k 



sin fc^ fjU^ca- 

Let us focus in the Dyson's equation with coordinates 



lea, I' , belonging to the central system: 

n 

I" 

n 



(A4) 



where I" runs over all the spacial indexes of the cen- 
tral system while kn,a labels degrees of freedom of the 
reservoir represented by Ha- The Green's functions with 
mixed coordinates kn,a, I' , in turn, satisfies the following 



equation: 

-Wa,kGl^,,^,{uj) - AaF^^ ^^i,^,{oj) = 0, 

+Wa,kFil^,,A^) ~ A*aGl ^^i,Ju;) = 0. (A5) 

After some algebra, the above equations can be casted as 
follows: 

-Wa,kFil^,,,A^)], (A6) 



+F;^'^iL.)F/l^,^Au;)i (A7) 



with: 



9KJiA 



uJ + efe„_„ + if] 

{u) + ekn,a) 
Aa 

with 77 = 0+ and E'^{ek^ J = 4„ „ + ^a- 



(A8) 



Substituting (jA6p into (jA4p . the latter equations can 
be expressed in the following way: 

-Y.^i...,i"G?',vAA^K..,i', 

I" 

[u: ~ Y.fJ'^\u:)]Fl^_,yLo) + E/«-«(c.)G«^,, ,(^) 

+ E^'-.'"^A;',c.H = 0- (A9) 
I" 

Notice that all the spacial indexes of the above equations 
run over coordinates of the central system, while the in- 
dexes corresponding to the reservoirs have been elimi- 
nated by defining the 'self-energies': 



E 



Wa.k\ 



, (AlO) 



being X"^" {u),ek„A = ^i^.i^'i^ ±£fc„,„) for ly = 5,/, re- 
spectively and A^'-'' (u;,£fc„.„) = ['^'''''^(w, £fc„_„)]* = A^. 

These steps can be repeated with each contact, which 
allows for the one by one elimination of the degrees of 
freedom of all the wires. The limit to the size of the 
wires going to infinite is summarized in Appendix (B) 



APPENDIX B: RETARDED SELF-ENERGIES 
ASSOCIATED A ID S-WIRE. 

We now evaluate the spectral functions Ta"^ (uj) = 
—2lm[Ila'^ '^{A] corresponding to the self-energies de- 
fined in the previous appendix in the thermodynamic 
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limit, Net — * oo. This corresponds to replacing 
(Na/Tr) dk in the expressions (jAlOp : 



(^) 



2 (•2u)q— /J 



E{u) 

x{5{uj-E{u))-5{ijj + E{u))}. (Bl) 



The final result is: 



x{[tj + r(w)]s+(w) + [tj - r{Lo)]s~ {uj)} 



x[s+(w) + s-(w)], 



2w'^ r{uj) 



(B2) 



with r(w) = e(|tj| - |Aa|)^w2 _ ^2 j^^^j ^ij-^) ^ 
0(|2wq| — \r{Lj) ± ^|)^4u'^ — (r(cj) ± It can be ver- 
ified that, for /i = 0, r^^(a;) reduces to the Im of the 
diagonal component of the self-energy defined by an infi- 



nite tight-binding wire with local pairing reported in Ref. 

[H. 

The final expressions for the retarded self-energies in 
the thermodynamic limit can be obtained by recourse to 
the Kramers-Kronig relation: 



+°° duj' T''^'''{lu') 
2Tr UJ — uj' + irj 



(B3) 



APPENDIX C: DYSON'S EQUATION FOR G< 
AND F<. 



The lesser counterpart of ([8]) is: 

[lu; - t^^'^iLo) e]G<(c.) + S^^^«HF<(c.) 

= tff'<iu)F^{Lo) - tf^^<iL0)G^{L0). (CI) 
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